Good Quantum Convolutional Error Correction Codes And Their Decoding 
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Quantum convolutional code was introduced recently as an alternative way to protect vital 
quantum information. To complete the analysis of quantum convolutional code, I report a way 
to decode certain quantum convolutional codes based on the classical Viterbi decoding algorithm. 
This decoding algorithm is optimal for a memoryless channel. I also report three simple criteria 
to test if decoding errors in a quantum convolutional code will terminate after a finite number of 
decoding steps whenever the Hilbert space dimension of each quantum register is a prime power. 
Finally, I show that certain quantum convolutional codes are in fact stabilizer codes. And hence, 
these quantum stabilizer convolutional codes have fault-tolerant implementations. 
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I. INTRODUCTION 



Quantum error correcting codes (QECCs) and their 
fault tolerant implementations are effective ways to pro- 
tect and to manipulate quantum information in the pres- 
ence of noise. A QECC works by adding suitable re- 
dundancy in the form of entanglement to the original 
quantum state in such a way that one can reconstruct 
the original state after decoherence. Since the discov- 
ery of QECC by Shor researchers have discovered 
many ways to construct QECCs. (See, for example, 
Refs. Ig-Q.) These QECCs share a common charac- 
teristic, namely, one first divides the original quantum 
registers into separate blocks of a fixed finite length. One 
then applies the same encoding method to quantum reg- 
isters in each block. Hence, this kind of codes are called 
quantum block codes (QBCs). A QECC can be decoded 
by first measuring the error syndromes of the encoded 
quantum state, and then by applying a necessary unitary 
transformation to the corresponding erroneous quantum 
registers For QBCs, this can be carried out in a 

block by block basis. Since there is only a finite number 
of error syndromes and hence also a finite number of re- 
covery operations in each block, decoding a QBC requires 
only a finite amount of work per (decoded) quantum reg- 
ister. 

Recently, Chau constructed another class of codes, 
known as quantum convolutional codes (QCCs), whose 
encoding scheme for each block depends on the states of 
several other blocks iQ. For example, he showed that 
the QCC 
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where ki G "Ln for all z > 0, /cj = for all j < 0, ujn is a 
primitive A'^th root of unity, the sum is from to A— 1, 
and all additions in the state ket are modulo A^, is capa- 
ble of correcting one error out of every eight consecutive 
quantum registers. 

While QCCs are of interest of their own right, it is not 
clear how to decode them. This is because the length 
of the original quantum state and hence the number of 
correctable errors by the code may both be infinite. Fur- 
thermore, decoding errors may propagate from one block 
to another due to their convolutional nature. Besides, it 
is not obvious how to manipulate a QCC fault-tolerantly. 

In this Paper, I address the questions of decoding and 
fault-tolerantly manipulating a QCC as well as a con- 
dition for the existence of a QCC whose decoding error 
does not propagate indefinitely. My key observation is 
that many classical convolutional code as well as quan- 
tum block code concepts can be extended to the quan- 
tum case when one performs the relevant operations care- 
fully. I first show that the well-known Viterbi decod- 
ing algorithm (VDA) for classical convolutional 
codes can be generalized to QCCs. Then, I show that 
the quantum version of Viterbi decoding (QVA) is equiv- 
alent to the maximum likelihood decoding. And hence, 
QVA is optimal in a memoryless channel. After that, 
I investigate the decoding error propagation in QCCs. 
In particular, I prove three equivalent criteria for QCCs 
to have finite decoding error propagation whenever the 
Hilbert space dimension of each quantum register is a 
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prime power. And finally, I address the question of fault- 
tolerant manipulation of QCCs by showing that the well- 
known fault-tolerant stabilizer code theory can be gener- 
alized to QCC. 



II. CLASSICAL AND QUANTUM VITERBI 
DECODING ALGORITHMS 

Before I go on, let me stress that in my subsequent 
discussion, I shall restrict myself to consider those classi- 
cal or quantum convolutional codes whose encoding can 
be implemented by a /c-input n-output (and hence also 
n — k preset registers in the quantum case) m-memory 
(that is, the encoding scheme depends on the state of the 
previous m blocks) quantum sequential circuit. (Com- 
pare with the definition of classical convolutional codes 
in Ref. ||l^.) Actually, all useful convolutional codes be- 
long to this category. 

Now, let me begin by briefly reviewing VDA for binary 
signals [p^-[l7|] . The algorithm starts by computing the 
Hamming distances between the first dn bits of the signal 
in the encoded sequence with the first dn bits of the 2^"^ 

possible code branches where d = -\- ij . One then 

keeps only those 2^'^"^^'^ code branches with small Ham- 
ming distances. (In case of a tie in Hamming distances, 
one keeps the corresponding code branches arbitrarily.) 
Now, one computes the Hamming distances between the 
first {d+l)n bits of signal with the first {d+l)n bits of 
all possible code branches that are consistent with a pre- 
viously kept code branch. This process is repeated until 
either the signal ends for signals of finite length or the 
process is repeated definitely for an infinitely long signal. 
The final surviving code branch is the decoded signal. In 
essence, VDA tries to find a codeword with the smallest 
Hamming distance from the signal ||l5|-pl[|. 

Clearly, VDA in the above form cannot be applied di- 
rectly to QCCs as it requires a complete knowledge of the 
encoded signal. However, if one examines the algorithm 
carefully, it becomes clear that what really required are 
the error syndromes of the encoded signal for different 
possible code branches. Consequently, VDA can be ap- 
plied to quantum signals. Suppose TZ spans of the set of 
all error recovery (unitary) operators for every n consec- 
utive quantum registers. (For a general error correcting 
code, TZ can be chosen to be in the form (^j^j^ (/i o Si) 
where fi is either an identity or a controlled phase shift 
operator on a single quantum register, and si is either an 
identity or a spin permutation operator on a single quan- 
tum register. Consequently, there are {N Nl)'^ elements 
in TZ.) Recall that error syndromes can be regarded as 
operators whose actions have no effect on an error-free 
encoded quantum state. (For instance, the action of er- 
ror syndrome on a stabilizer code simply permutes the 
stabilizer.) Thus, by measuring the eigenvalue of an er- 
ror syndrome, one gains some information on the location 
and kind of error occurred in a quantum code |^,^ . Since 



the set of all error syndromes is closed under composition 
and that the QCC is of finite memory, there are only a 
finite number of independent error syndrome operators 
that acts only on the first d blocks of encoded quantum 
registers. So, for each R E TZ^'^, I use a finite (and fixed) 
number of ancillary quantum registers to measure the 
error syndromes of the first dn quantum registers after 
subjecting to the unitary operation R. In this way, I can 
the erroneous registers for each R e 7?,®''. 

Once the erroneous quantum registers are located, how 
are we going to correct the quantum errors? To answer 
this question, I have to introduce the following definition 
first. 

Definition 1 Let {"^i) and {"^2) be two quantum signals 
of possibly infinite lengths. If it is not possible to find 
a unitary operator involving a finite number of quantum 
registers which maps j^E'i) to |5'2), then I say that the 
quantum Hamming distance between and |\E'2) 
is infinite. Otherwise, I define the quantum Hamming 
distance between these two quantum signals as the min- 
imum number of quantum registers involved in unitarily 
transforming from one state to the other. 

Similarly, I define the minimum quantum Hamming 
distance between a quantum signal \^) and the set of 
all possible codewords of a QECC C to be infinite if the 
quantum Hamming distances between \^) and all code- 
words of C are infinite. Otherwise, I define the quantum 
Hamming distance between \ and C to be the minimum 
possible quantum Hamming distance between \ ^) and the 
codewords in C . And for simplicity, I shall simply call 
the minimum quantum Hamming distance between \^) 
and C the quantum Hamming distance 0/ j^*). 

/ also define the recovery cost of bringing the first 
dn quantum signals with respected to R E TZ'^'^ to be the 
quantum Hamming distance of the quantum signal plus 
the (minimum) number of registers affected by R. 

Readers can easily check that quantum Hamming dis- 
tance between two quantum signals is a metric for the set 
of all quantum signals. Moreover, in a loose sense, the 
recovery cost measures how close and how much work is 
required to bring a quantum signal to a quantum code- 
word. 

With the above definition, I am ready to report QVA 
decoding: By carefully measuring the error syndromes 
using ancillary quantum registers, I compute the recovery 
cost for the first dn quantum registers for each R S TZ'^'^. 
I keep |7?.|'^~^ error recovery operators with small recov- 
ery costs out of the \TZ\'^ possible ones (where \TZ\ de- 
notes the number of elements in the set TZ). Then, I 
go on to compute the recovery cost for the first {d-\-l)n 
quantum signals with respected to the set of all possible 
recovery operators in 7^®(''+i) that are consistent with a 
previously kept recovery operator in TZ'^'^. For a quan- 
tum signal of finite length, I repeat this process until 
the quantum signal terminates. Then, I regard the error 
of the signal to be caused by the one that produces the 
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minimum possible recovery cost among those \R\'^ ones I 
kept at the end. And I correct the quantum signal accord- 
ingly. For quantum signals of infinite length, I have to 
repeat the recovery cost selection process infinitely many 
times in order to find the minimum recovery cost path. 
In practice, we usually run the QVA over a large but fi- 
nite number of quantum blocks and decode the signal in 
each block separately. The length of such quantum block 
is usually adaptive; that is to say, it is chosen in such a 
way that the recovery cost paths retained do not differ 
very much from each other. In this way, the effect due 
to the choice of the length of the block is minimized. In 
summary, regardless of the length of the quantum signal, 
QVA uses a finite number of operations on average to 
recover it. 

Finally, I need to convert the recovered encoded signal 
to its unencoded form. Since I may have an infinitely long 
signal, the usual trick of running the reversible encoding 
quantum circuit backward does not work. Hence, I have 
no choice but to decode the signal starting from the first 
encoded block. Remember that by including the preset 
quantum registers, the encoding process can be repre- 
sented by a unitary transformation. Let C is a fc-input 
n-output TO-memory QCC and be a quantum sig- 
nal (with preset quantum registers added). I denote the 
encoding process for this quantum code C and quantum 
signal 1^) by Ui. Moreover, I denote the encoding pro- 
cess by the same code C on the quantum signal by Ui 
where equals l^*!) except that the first k{i~l) quan- 
tum registers are set to zero. Using this notation, I write 
the encoding process 

Consequently, the decoding process is given by 

f/f' = • • • o (c/g o u^^y^o {U2 o uyy^o {Ui o uyy' 

= • • • o (C/4 o o ({/3 o {/2-I) o (c/2 o C/-1) . (2) 

Note that for the m-memory code C, Ui+i o U^^ is a 
unitary operator acting only on the ith, (i-l-l)th, up to 
(j-l-m)th encoded blocks for all i. Moreover, it is easy to 
check that the action of U2 o U^^ is to extract the state 
of the first unencoded block of quantum registers out of 
the encoded state. After that, the action of U3 o U^^ is 
to extract the state of the second unencoded block out 
of the encoded state, and so on. Thus, Eq. (|^) gives us a 
way to decode the QCC C. 

Example 1 It is easy to compute Ui+ioU^^ in practice. 
For example, to decode the QCC in Eq. (yj in Eq. 
I subtracts the second quantum reqister by the first, the 
fourth by the first and the third, the fifth by the first, the 
sixth by the first and the fourth, the seventh by the third, 
and finally the eiqhth by the third. The resultant quan- 
tum state is qiven by Eq. (^^. Then, I discrete inverse 
Fourier transform the first quantum register. The resul- 
tant state is given by Eq. after summing over the 
dummy index pi . Finally, it is straight forward to uni- 
tarily convert the state in Eq. to the state in Eq. 
by multiplying a phase proportional to the product of the 



first quantum register by the sum of the third and the 
twelve registers, and then followed by discrete Fourier 
transforming the third quantum register. The result of 
all these steps above gives U2 o Ui^. 
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In this way, I decode the first quantum register from the 
code using a finite number of two body operators. And 
inductively, I can decode the rest of the encoded quantum 
signal efficiently. 

Now, I move on to prove the optimality of QVA for a 
memoryless quantum channel, namely, a noisy quantum 
channel whose error occurs randomly and independently. 
Similar to VDA, QVA tries to search for a solution in the 
codeword space with a minimum recovery cost from the 
signal. (I choose to minimize the recovery cost instead of 
simply the quantum Hamming distance because certain 
fault tolerant operations U for the QCC may belong to 
0+J;T ^- In this case, the quantum Hamming distances 
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of \^) and U\^) agree although they have different re- 
covery costs.) After transmitting an arbitrary unknown 
encoded quantum state through a memoryless channel, 
the effect of decoherence can be regarded as a Markov 
process. In other words, the probability that the error 
recovery operator needed to act on the (t + l)th block of 
quantum signals, given all error recovery operators for 
the first t blocks, depends only on that of the (t-l-l)th 
block. Since I do not manipulate my encoded quantum 
signal during its transmission, I can always conceptually 
regard my error syndrome measurement to be performed 
immediately after the errors are introduced into the sig- 
nal. But once I have measured the error syndromes, the 
location as well as the type of error each quantum register 
is suffering becomes classical data. Therefore, the effect 
of a quantum memoryless channel is the same as that 
in a classical probabilistic memoryless channel. More 
precisely, I can always model the chance for a certain 
quantum error R G TZ to error in a quantum register by 
a classical probability function. (Compare to the argu- 
ment used in the proof of the security of quantum key 
distribution in Re f. p2| .) Consequently, the optimality 
proof of VDA |15[-|l9| carries over directly to QVA. 

Similar to classical convolutional codes, there are two 
important probabilities which measure the performance 
of a QCC. The first one is called the error probabil- 
ity Pe{E), which is defined to be the probability that 
a wrong decoding path is chosen at any given timestep. 
And the second one is called the qubit error probability 
Pfc (E) , which is defined to be the expected number of in- 
formation qubit decoding errors per decoded information 
qubit. For A'^ = 2 in a binary symmetric channel, these 
two probabilities are given by p7t] 



Pe{E) < Ad2V 



(4a) 



becomes exponentially small. Thus, we have a family of 
good QCCs. 



III. NON-CATASTROPHIC QUANTUM CODES 

The ability to decode a QCC is not sufficient to make 
QCC useful. We must also make sure that any decoding 
error will not propagate infinitely in spite of the convo- 
lutional nature of the code. To facilitate discussions, I 
borrow the following terminology from classical coding 
theory: 

Definition 2 A QCC is said to be catastrophic if there 
exists a local decoding error that can propagate infinitely. 
Otherwise, a QCC is said to be non-catastrophic. 

Clearly, useful QCCs must be non- catastrophic. 

In case of classical convolutional codes and when the 
number of internal states per register A'^ is a prime power, 
a convolutional encoder can be mathematically repre- 
sented by a polynomial of one variable over a finite field. 
Such a polynomial ring is clearly an Euclidean domain. 
In particular, two polynomials in an Euclidean domain 
have a unique greater common divisor (up to multiplica- 
tion of units). Using this nice property of an Euclidean 
domain, Massey and Sain |25 2^ proved a necessary and 
sufficient condition for a classical convolutional code to 
be non-catastrophic. Nonetheless, quantum mechanical 
operations are intrinsically non-commutative. Thus, the 
proof of Massey and Sain does not work for QCCs. 

Quite surprisingly, a necessary and sufficient condition 
for a QCC to be non-catastrophic can still be found when- 
ever N is a prime power (and hence Z^r is a finite field). 
And I am going to report the criterion after introducing 
the following rather involved notations. 
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where d is the minimum quantum Hamming distance be- 
tween codewords. Ad is the number of mutually orthogo- 
nal encoded states of quantum Hamming weight d, Bd is 
the number of non-zero mutually orthogonal information 
qubits on all weight d paths, k is number of information 
qubits per block and p is error probability of the channel. 

For a fixed QCC, the distance of the code d is finite 
and hence both Pe and Pb scales only like a power law of 
p. Nevertheless, there are many fc-input n-output QCCs 
with different memories m. And for many fixed k and 
n, the distance of the code d increases approximately lin- 
early with memory m. And such a family of codes may be 
constructed, for instance, from a corresponding family of 
classical convolutional codes. One such family of classical 
code as briefly discussed in Ref. |0 encodes one classical 
bit into two. By the construction of Chau in Ref. it 
can be turned into a family of 1-input 4-output QCCs. 
And as the memory m tends to infinity, both Pg and Pb 



A. Notation For Encoding And Decoding Qubits 
When A Is A Prime Power 

In case A'^ is a prime power, any two body unitary 
operations can be generated by the span of the follow- 
ing elementary Pauli group operations p^ : (a) addition 
|a;) ^ \x + a) and \x,y) i~+ \x + y,y) for some a £ Zin, 
(b) multiplication |a;) i— > \ax) for some a G ZAr\{0}, (c) 
Fourier transform \x) i— > X^y'^jv'ly)' i'^) local phase mul- 
tiplication \x) ^ Lo'ff\x), and (e) non-local phase multi- 
plication \x,y) ^ uj'j^^\x,y). 

Let me first considers those QCCs that can be en- 
coded can be implemented by a fc-input n-output m- 
memory quantum sequential circuit. In this case, I can 
group the initial unencoded quantum state and preset 
registers into blocks of length n. This state is spanned 
by {(8)£T N»i'2^i2, ■ ■ ■ ,3;™)}, where for i > 1 and 
^ < i < n — k are the preset registers and Xij for i > I 
and n— k < j < n are the quantum information reg- 
isters. With this notation in mind, I define the follow- 
ing operators: (a) state delay operator D"^ : \xij) i— > 
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), (b) Fourier transform operator Fp 



^ijViJ — p 

N 



Hij-p), (c) phase projection operator^ P 



(d) local phase multiplication operator 
L™ : \xij) I— > u/^^'^\xij)^ (e) non-local phase multiplica- 
tion operator Mjnp '■ i-^ a;^'''^'~"''^~''|xij-), (f) state 
addition operator + : {e''''\x),e'^'\x')) ^ e'('^+"*')|a; + x'), 
and (g) state multiplication operator 
Using these notations, the five elementary Pauli group 
operations can be represented by (a) l + D^P, (b) a, (c) 
Fq, (d) L", and (e) M^q, respectively. More generally, if 
I write the initial state ket together with the preset and 
ancillary quantum registers in apx 1 column vector, then 
the composition of several Pauli group operations can be 
represented hy a, p x p matrix whose elements belong 
to the non-commutative ring Zn {D, Mij, P, L, Fi) with 



F' 



4, = P, PD = DP, LM^, = MuL, PL = P, 



PM,j = P, AhjMpq 



MpqMij, and Mijm 



'Mr 



for m e Zn\{0}. 

Since the Pauli group spans the set of all two body op- 
erators |2^, a general quantum encoding circuit J/encodc 
for a fc-input n-output m-memory QCC can be writ- 
ten as a finite sum '^j_{cti,gi) where € C and 

e {Zn (P',Mjj,P,P,P,))"+^'"+^ where p is the num- 
ber of ancillary quantum registers required in the en- 
coding process per block. For instance, the operator 
E»=o^(V^, L') sends |0) to |0) and all other \i) to 0. Fur- 
thermore, the unitary operator sending \x,y) to \x,x+y) 
can be written as 



(1,1) (0,0) 
(1,^^) (1,1) 



(5) 



Readers should observe that the phase projection opera- 
tor P in Eq. (^ is essential. If I replace P by 1 in Eq. (||), 
then the replaced operator will not be well-defined for it 
would have mapped e'^^^~^^^\x,y) to e'^^^'^^^\x,x-\-y) and 
e'*i|x) e^^^ly) to e't^^i+^^lja;, x-l-y). In addition, the 
operator expressed in Eq. (|5) is unitary in spite of its 
apparent non-skew symmetric form. 

Let me denote the set of all finite sums in the form 
^j(ai,gi) by K. Then, if I forget about the initial 
preset and ancillary quantum registers and simply rep- 
resent the initial unencoded quantum information as a 
fc X 1 column vector, then I can simply write a fc-input n- 
output m-memory quantum encoding circuit as a n x fc 
matrix in K"'''. The decoding circuit for this QCC is 
equal to a {n+p) x (n+p) matrix U~^^^^^. Nevertheless, 
U, ' 



encode 



circuit, if I forget about the initial preset and ancillary 
quantum registers used in the decoding circuit, then I 
can present the decoding circuit by a fc x n matrix. 

Example 2 Using the above notation, the encoding and 



decoding algorithms for the classical non- catastrophic 
convolutional code (written in a quantum state ket form) 



|fci,/c2, • ■ •) 
written as 



and 



\ki + ki_2,h + ki-i + ki-2) 



{1,P[1 + D^]) 
{1, D + P[l + D^]) 



be 



(6) 



(7) 



respectively. 

Similarly, the encoding and decoding algorithms for the 
classical catastrophic convolutional code \ki, k2, . . .) ' — > 
'S'iJo \ki + ki-i, ki-\-ki^2) can be written as 



{1,P[1 + D]) 
{1,1 + PD^) 



+ PD]-^D) + PD]-^)] , 

respectively. 



(8) 



(9) 



Example 3 One possible way to encode the quantum 
state \ki, k2, . . .) as a QCC given by Eq. ^ is as follows: 
First, I prepare a number of preset quantum registers and 
write the initial state as (^^tjT I'*'*' ^' ^' ^) ■ Then, I trans- 
form this state to {^^t!T l^i: ^i-i? 0) by the unitary 
operator 



(1,1) (0,0) (0,0) (0,0) 

{l,PD) (0,0) (0,0) (0,0) 

{1,PD^) (0,0) (1,1) (0,0) 

(0,0) (0,0) (0,0) (1,1) 



(10) 



Then, I transform the state to ^^J^ \ki+ki^2, ki+ki^i- 
fc,_2,0,0) by 



(1,1) (0,0) {l,PD) (0,0) 

{l,PD) (1,1) {l,PD^) (0,0) 

(0,0) {l,-PD) (1,1) (0,0) 

(0,0) (0,0) (0,0) (1,1) 



(11) 



^ i=l 



Next, I unitarily transform 

+ 00 pi{ki+ki-2)+gi{ki+ki-i+k, 

'■'jv 



the state 
\pi,q.i, 0,0) by 



to 



(l,Po) (0,0) (0,0) (0,0) 

(0,0) (l,Pi) (0,0) (0,0) 

(0,0) (0,0) (1,1) (0,0) 

(0,0) (0,0) (0,0) (1,1) 



(12) 



shall explain why I introduce such a non-invertible pro- 
jection operator in the next paragraph. 
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Finally, I bring the state to Eq. hy the unitary trans- 
formation 



{l,P[l + D]) (0,0) (1,1) (0,0) 

{l,P[l + D]) {1,PD) (1,P) (1,1) 

(0,0) (1,1) (0,0) (1,P) 

(1,1) (1,P) (0,0) (1,P) 



(13) 



Thus, the unitary encoding transformation J/oncodc for the 
QCC in Eq. simply equals to the product of the ma- 
trices in Eqs. (|7^-([7^. Moreover, if we forget about 
the initial preset registers, then the encoding operation is 
simply given by the first column of the matrix J/encode, 
which is given by 



{l,P[l + D]Fo[l + PD^]) 
(1, P[l + D]Fo[l + PD^] + PDFi [1 + D + D^]) 

+ PD]FiP[l + D + D^]) 
(1, Po[l + PD^] + P[l + D]FiP[l + D + D^]) 



(14) 



Similarly, the decoding operation is equal to the first row 



of the matrix U 



ode' 



namely, 



{l,P[DF,-'-l]) 
{1,-PDFi) 
(1, [1 + PD]FoP - PDF^^ - P[l + D]) 
{l,[l + PD]F-'+P[l + D]) 

{1,P[~DF,{-1)~1]) 
il.-PDFi) 
(1, [1 + PD]FqP + PDPi(-l) - P[l + D]) 
{l,-[l + PD]Fo{-l) + P[l + D]) 



(15) 



B. Criterion For Non-Catastrophic Quantum Code 
When TV Is A Prime Power 

Now, let me report a useful lemma before proving a 
necessary and sufficient condition for non-catastrophic 
QCCs. 

Lemma 1 Suppose N is a prime power. And let M G 
K^'P be a valid unitary operator acting on a possible in- 
finitely long quantum signal. Then, M can be decomposed 
into a product of finite product J^^^j^ Mi. Moreover, the 
P' 

for all i 



^ elements in each matrix M; commute with each other 



Proof: Since is a prime power and hence Zjv is a field, I 
can borrow the idea in Ref. [p4[ to decompose the matrix 
M as a product of finitely many matrices. Observe that 
I can always find an invertible p x p matrix A'^i such that 
the element located in the pth row and (p— l)th column 
in N^^M equals (0, 0). Besides, I can choose iVi in such 
a way that the elements in the ith row and jth column 
satisfy the condition 



(iVi 



(1,1) if i = j and i < p— 2 

(0,0) if and (i,j) ^ (p,p~l) 



(16) 



Similarly, I can find a p x p matrix 7V2 such that 
(N^'^N^'^M).. = (0,0) whenever = {p,p-l) and 

{p,p~2). Moreover, 



iN2 



(1, 1) if i = j and i 7^ 2 or p 
(0,0) if i and (z,j) (p,p-2) 



(17) 



Inductively, I 



N 



p{p-i)/2y[p{p- 
lar matrix. 

N, 



can find Ni such that M' = 
i)/2]-i ' ' ' ^r^^^ is a upper triangu- 
Besides, at most three elements in 



Ip are non-zero where Ip denotes the iden- 
tity operator. Similarly, I can transform the ma- 
trix M' into a diagonal one by means of p{p—l)/2 
matrices in a similar form as Ni. Thus, M = 

N1N2 ■ ■ ■ Np(^p_iy2M"Nipj^p_iy2] + lN[p(p-i)/2]+2 ' ' ' ^p(p-l) 

where M" is a diagonal matrix and all Ni — Ip can be 
brought into the form in Eq. ( p^ by relabeling some 
columns and rows plus possibly a transposition. 

It is obvious that M" is equal to a product of p diagonal 
matrices, each of which has at most one diagonal element 
differ from (1,1). Besides, elements in each of the p ma- 
trices commute with each other. Thus, to complete the 
proof, it remains to show that Ni can be decomposed into 
a finite product of matrices whose elements commute. In 
fact, it suffices for me to decompose it for the 2x2 matrix 
\ A B 1 

M = (0 0) C ' decomposition for the matrix Ni 

is similar. Since M is a well-defined operator, either A 
or B (but not both) must be in the form PX for some 
X G K. In the first case, A = PX and it is easy to check 
that 



M = 



(1,0) 
(0,0) 



(0,0) 
C 



(0,0) 



B 
(1,0) 



X 
(0,0) 



(0,0) 

(1,1) 



(18a) 



And in the second case, B = PY and 



M 



(1,0) (0,0) 
(0,0) C 



(1,P) Y 
(0,0) (1,0) 



X (0,0) 
(0,0) (1,1) 



(18b) 



Since element s in each o f the matrices in the right hand 
sides of Eqs. (18a) and (18b) commute, so the lemma is 
proved. □ 

After going through the above preparatory discussions 
and examples, I am ready to state a necessary and suffi- 
cient condition for non-catastrophic QCCs. In fact. The- 
orem |l| below generalizes a necessary and sufficient con- 
dition for classical non-catastrophic codes pqjpq] . 
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Theorem 1 The following statements concerning a It- 
input n-output m-memory QCC are equivalent when N 
is a prime power: 

(a) The QCC is non-catastrophic. 

(b) There exits a quantum encoding circuit (that in- 
cludes the preset and ancillary quantum registers) g G 
j^n+p,n+p ^^jg Q(jQ gy^^j ^f^^^ 1^1^ inverse g'^ 

exists. Moreover, elements in the matrix g~^ can 
be expressed as a finite sum '^i{<^i,g'i) with g[ € 
In {D, D'-^,Mpg, P, L, Ep) for all i. 

(c) There exists a quantum encoding circuit that can 
be decomposed into the following finite product g — gi 
in such a way that for each i, (1) gi € j^n+p,n+p^ 
elements of matrix gi belong to a commutative poly- 
nomial ring, and (3) the inverse (detgi) ^ exists and 
can be expressed as a finite sum '^i{oij,g'j) with g'j € 
ZN{D,D-\Mpg,P,L,Fp). 

(d) The quantum encoding circuit (that excludes the 
preset registers) h S K'^'-'^ can be expressed as a finite 
product of matrices hi in such a way that each i, 
(1) hi G ^ (2) elements of matrix hi belong to a 
commutative polynomial ring, and (3) the greatest com- 
mon divisor ti of the determinant of all the sub- 
matrices of hi is invertible and t^^ = with 
t'^^ eZN{D,D-\Mpq,P,L,Fp). 

Proof: By suitably adding ancillary quantum registers as 
well as enlarging the encoding matrix to include those 
ancillary registers, it is easy to see that (d) (c) =J> 
(b). Now, I move on to show that (b) (a). Since 
there is a decoding circuit that can be represented as a 
{n -\- p) X (n +p) matrix whose elements X^iC'^ii ^i) with 
h, e In {D, D-\M^j,P, L, F,). In other words, decod- 
ing each quantum register in the code requires only in- 
formation from a finite number of encoded quantum reg- 
isters. Thus, if there are only a finite number of encoded 
quantum registers in error, then the decoding errors will 
only be localized in a finite number of quantum registers. 
Hence, the code is non-catastrophic. 

To complete the proof, it remains for me to show that 
(a) ^ (d). Recall that if ?7oncodc S iC"'" is the encoding 
circuit, then the decoding circuit equals C4~codo- 
statement (d) is false, then I can extend the fc x n decod- 
ing circuit into a n x n one. And since iV is a prime power, 
so by Lemma 0, I can conclude that elements in the fc x n 
decoding circuit are in the form ^i) where hi be- 

longs to the formal power series non-commutative ring 
Zn {{D, Mij, P, L, Fi)) but not every element in the de- 
coding circuit belongs to Zn {D, Mij , P, L, Fi) . Conse- 
quently, there exists an encoded quantum register whose 
state aff'ects the states of infinitely many decoded quan- 
tum registers. Thus, the QCC is catastrophic and this 
complete the proof. □ 

Now, it is clear from the proof of Theorem Q that if 
I first let the encoded quantum to go through the QVA 
and then I apply the unitary transformation g~^ (which 



is the left inverse of 5) to it, I can recover the original 
unencoded quantum information. In addition, it is also 
clear that any QCC that cannot be expressed as a fc-input 
n-output m-memory sequential quantum circuit must be 
catastrophic. Moreover, the conclusion in Theorem ^ re- 
mains valid if I extend the meaning of m-memory QCC to 
include those QCCs whose encoding scheme depends on 
the state of a finite number of previous or future blocks. 

One possible way to construc t Q CC is to start with 
a classical convolutional code C ||l^ . Chau showed that 
one can first encode a quantum signal using the classi- 
cal code C, then one takes the local Fourier transform 
on each encoded quantum register, and finally one en- 
codes the resultant state ket by the code C again, one 
gets a QCC Here, I show that the QCC generated 
this way inherits the error propagation behavior from its 
parent classical code. 

Corollary 1 Suppose C is a k-input n-output classical 
convolutional code and Q be the corresponding -input 
-output QCC obtained using the above method. Then 
C is catastrophic if and only if Q is catastrophic. 

Proof: I write the quantum encoding scheme as a product 
of three matrices (715233 where gi and (73 involve the sym- 
bols D and DP and g2 involves the symbol Fi . That is to 
say, gi and 53 represent the initial and final encoding by 
the classical code C and 172 represents the local Fourier 
transform. Suppose C is non-catastrophic, then clearly 
I can arrange gi, 52 and 53 to satisfy statement (d) in 
Theorem |] p5| , p6t . Hence, Q is non-catastrophic. Con- 
versely, if C is a catastrophic code, then from the con- 
struction of Q, it is clear that one can always find a finite 
number of spin flip errors for Q such that the decoding 
errors propagate infinitely. Hence Q is catastrophic. □ 

Corollary implies that the QCC given in Eq. (0) is 
non-catastrophic. 



IV. FAULT-TOLERANT COMPUTATION USING 
QUANTUM CONVOLUTIONAL CODES 

The ability to decode a non-catastrophic QCC is still 
not enough to make them truly useful. We have to impose 
the requirement that the QCC must have a fault-tolerant 
implementation so that quantum information processing 
can take place in the encoded form. In QBCs, we know 
that all stabilizer (block) codes have a fault-tolerant im- 
plementations 27 -3^] under suitable wiring of quantum 
gates. And now, I am going to generalize the theory of 
stabilizer code and its fault-tolerant implementations to 
the world of QCCs. 

For stabilizer codes, I restrict myself to consider the 
case when N = 2. Recall that in the case of QBC and 
when = 2, if we denote the coding space of an n qubit 
code by T, then the stabilizer of this code S is some 
Abelian subgroup of the group 7?,®" whose elements fixes 
T P |^, |3l| , |3^ . Besides, S can be generated by a finite 
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number of operations gi € 7?.®" , known as the generators 
of S. FinaUy, the encoded spin flip and phase change 
operations are specified in 7?.®". These operations com- 
mute with the stabilizer S. More precisely, the codeword 
for a /c-input rt-output stabilizer QBC can be written (up 
to an overall normalization constant) as: 



\Xi,X2 



I ■^Zj ■ • • 7 



1 

E 



where gi and the generators of the stabilizer and 

the encoded spin flip operation for \xj), respectively 
In addition, the encoded phase shift operators ct^ j exist 
in 7^®" for aU j. 

Generalizing the stabilizer (block) code formalism to 
the QCC world is easy. One only needs to be more care- 
ful in dealing with the infinite number of qubits and hence 
the infinite number of generators for the stabilizer. First, 
one replaces TZ^"- by Clearly, ^^"^ hence 

the stabilizer S have a counterable number of generators. 
Thus, Eq. (|l|) holds for QCCs as A; ^ oo. Besides, for 
a rn-memory QCC, the encoded spin flip operators 
as well as the encoded phase shift operators ct^j act on 
no more than 0(n(m + 1)) qubits. In this way, the fault- 
tolerant error syndrome measurement procedure in sta- 
bilizer block code [^,0 directly applies to convolutional 
code. Finally, one concatenate the QCC with another 
stabilizer QBC to L levels. Then, by correcting the er- 
rors in all levels concurrently, one achieves an error re- 
duction from 0(e) to O(e^). Hence, Gottesman's ||3|] 
proof that all stabilizer codes have fault-tolerant imple- 
mentation directly carries over to the QCC world. (See 
also Ref. |^3j for related results.) Thus, non-catastrophic 
stabilizer QCC are good codes. Recently, Gottesman ex- 
tended his theory to cover a large number of A^ary fault- 
tolerant quantum codes using Pauli group p^ ]. A direct 
consequence of his result is that we can easily construct 
many A^ary fault-tolerant QCCs. 

Finally, I go on the show that the QCC in Eq. (|l|) is a 
stabilizer code. In fact, I prove something more general: 



Theorem 2 Let C be a classical convolutional code. And 
let Q be the corresponding QCC as described in Corol- 
lary y. Then, Q is a stabilizer code. 

Proof: The proof follows directly from the three lemmas 
below. □ 

Lemma 2 All classical binary (block or convolutional) 
codes are stabilizer codes. 

Proof: Without lost of generality, I consider a m-memory 
classical convolutional code. Then, the encoded spin flip 
operator Wxj is nothing but a finite number of ax act- 
ing the encoded qubits. Since the code is classical, the 



encoded state for each \xi,X2,. . .) can simply be repre- 
sented by a single state ket without any dummy summa- 
tion index. More precisely, elements of the stabilizer are 
those commute with the encoded spin flip operators and 
at the same time can be expressed in the form AiA2At, ■ ■ ■ 
where Ai acts on the ith encoded qubit and Ai G {1, cr^}. 
Clearly, this kind of elements forms an Abclian subgroup 
of JlitiT ^ ^ counterable number of generators. 

Hence, the lemma is proved. □ 

Example 4 The stabilizer associated with the classical 



(19) block code \k) 



\kkk) is generated by UyUy^ and 



(Jylay 



Furthermore, the encoded spin flip and phase 
shift operators equal ax<Jx<^x and a z , respectively. 

Example 5 The encoded spin flip operators for the 
classical convolutional code |fci,fc2,...) i — > + 
ki_2,ki + h-i+ki^2) is given by axCfxtaxO-xCFxtt ■ ■ ■ , 

and so on. The encoded phase shift opera- 

tors are given by UcTzO'zU---, TLTLHILaz'^z^^ ' ' 
and so on. Besides, the stabilizer for this code 
is generated by cr^CTzll---, IcTzltTxCTzCrzll • • 
azO-ztaztazCTzO-zU ■ ■ ■, tta zU zta zta zO z(t z'^'^ ■ ■ 
IIIIct^CTzIctzIctzCTzCTzII • • and so on. 



Lemma 3 Let C be a classical binary code. And let C 
be the code obtained by locally Fourier transforming each 
qubits in the code C. Then, both C and C are stabilizer 
codes. 

Proof: Lemma ^ says that C is a stabilizer code. Suppose 
gi are the generators of the stabilizer of C, Wxj is the en- 
coded spin flip operators of C as described in Lemma ||. 
Define g'^ to be gi with ax replaced by <7z ■ Similarly, I de- 
fine axj to be a'x j but with ax replaced by az- Then, it 
is easy to verify that g'i generate the stabilizer of C". Be- 
sides, a'x j is the encoded spin flip operator for the code 
C. ' □ 



Lemma 4 Let C and C be the codes 



described 



Lemmas ^ and |^. Then the code C" obtained by first 
encode the state by C and then by C is a stabilizer code. 
( Compare with Ref. JS^/ for a similar result.) In addi- 
tion, each encoded spin flip and phase shift operator for 
C" acts on a finite number of qubits provided that C and 
hence C" are non- catastrophic. 

Proof: Suppose C and hence also C are fc-input n-output 
codes with finite memory. Then from Lemmas ^ and 
I can write the generators of the stabilizer code C as 
AiiAi2Ai3--- where Aij G {l,ax}. Moreover, I write 
BiiBi2 • • • as the generators of the stabilizer code C" 
where Bij S {IjCTx}. Suppose XiiXi2 ■ ■ ■ and Zi\Zi2 ■ ■ ■ 
be the encoded spin flip operators for codes C and C", 
resp ectively. Recall that C can be expressed in the form 
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+ 00 



\xi,X-2,---) 



(20) 



Since C is of finite memory, the sum in each of the en- 
code qubits in Eq. ( pO| ) above is finite. More precisely, 
aij = or 1 and for each fixed i, only a finite number 
of Qij equals one. Consequently, the QCC C" can be 
expressed in the form 



\X1,X2, ■ . .) 



1=1 



/ . - N 

LPl,P2,-- 



(21) 



where 6^ = or 1 and for each fixed i, only a finite 
number of bij equals one. 

If C is catastrophic, its decoding circuit can be ex- 
pressed as a formal power series matrix. While if C is 
non-catastrophic, its decoding circuit can be expressed as 
a polynomial matrix ||2^. (See also Theorem]^.) Conse- 
quently, there exist aj £ {0, 1} such that a^'^a^^^a^'^ ■ ■ ■ 
is an operator acting on the codeword of C" whose result 
is to map Pi to Pi + 1 for all i. Besides, a'^'^'^ al'^^ a"^'^ • • • is 
an operator acting on the codeword of C" whose result 
is to multiply the codeword by a phase (—I)''' for all i. 
Similarly, there exists dij S {0, 1} such that a^'^a^'^ ■ ■ ■ 
is an operator acting on the codeword of C whose result 



IS an operator 



is to map Xi to Xi+1, and that crf'^cr 
acting on the codeword of C whose result is to multiply 
the codeword by a phase (—1)^' for all i. Furthermore, 
for each fixed i, only a finite number of Cy and dij equals 
one if C is non-catastrophic. 

Once I know how to add one to Xi and pi as well as 
how to add phases (— 1)"^' and (—1)''* to the codewords 
of C and C" in the previous paragraph, I can use them to 
construct the encoded spin fiip and phase shift operators 
for the code C". They are given by 



IZjdijCji Y.jdijCj2 Y.jdijCj3 
(Jz Oz Oz 



(22a) 



dzOzOzOzOxOx'^'^'yzOzUzdz'^'^ ■ ■ ■, tlllazCTzazCrzaxO-xtl 

Oz^z^zOz^^ ■ ■ ■, and so on. In addition, the en- 
coded phase shift operators are Gx'JxCTx^Ox'^OxUx^^ • ■ 
11 1 1 cr j^CTj; CT j;l(Tj^l (7 jjCTx 11 • • cind so on. 

According to the proof of Lemma ^, fault-tolerant com- 
putation is possible for all QCCs constructed using the 
method in Theorem |^. More importantly, if one starts 
with a non-catastrophic classical convolutional code C, 
then the fault-tolerant spin flip, phase shift, and con- 
trolled swapping for the QCC C" constructed in Theo- 
rem H can all be done in finite number of quantum gates. 
In fact, as long as I carefully wire my quantum circuit 
to prevent the spreading of quantum errors throughout 
all the qubits (see Ref. for the tips of how to do 
this), I can perform fault-tolerant quantum computation 
on this kind of QCCs. Suppose I have a quantum signal 
\xi, X2, . . .) and if I follow the fault-tolerant computation 
wiring rule, I may even perform computation between the 
ith and jth encoded qubits in the above signal provided 
that their encoded spin flip and phase shift operators acts 
on distinct places in the encoded signal. 



V. DISCUSSIONS 

In summary, I have generalized the VDA to QCC and 
have shown the optimality of QVA for a memoryless 
channel. In addition, I reported a simple way to test 
if a QCC is non-catastrophic. The key observation for 
all these is that a lot of the classical coding concepts can 
be "quantized" provided that one performs the relevant 
operations with care. Finally, I show that certain QCCs 
can perform fault-tolerant quantum computation. Since 
classical convolutional codes may be regarded as stabi- 
lizer codes and good classical convolutional codes exists, 
therefore I conclude that good QCCs and their decoding 
algorithm exist. 



and 



(Jx 



Ox 



(22b) 



respectively. 

After identifying the encoded spin flip and phase shift 
operations in C", it remains for me to find the generators 
of the stabilizer of C". First, by direct checking, I know 
that the operator cr^'^cr^'^ • • • belongs to the stabilizer of 
C" for all i. Then, similar to the proof of Lemma |2| I con- 
sider operators in the form Ai\Ai2 ■ ■ ■ with G {1, Oz\ 
that commute with the encoded spin flip, encoded phase 



shift and ot^^ot 



Now, I choose a (counterable num- 



ber of) generators amongst them. Then, the union of 
these operators and cr^'icr^*^ • • • generate with stabilizer 
of the code C". □ 

Example 6 When N ^ 2, the encoded spin 
flip operators for the QCC in Eq. nhj are 
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